High-growth-rate magnetohydrodynamic instability in differentially rotating 

compressible flow 
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The transport of angular momentum in the outward direction is the fundamental requirement 
for accretion to proceed in an accretion disk. This objective can be achieved if the accretion flow 
is turbulent. Instabilities are one of the sources for the turbulence. We study a differentially 
rotating compressive flow in the presence of non vanishing radial and azimuthal magnetic field and 
demonstrate the occurrence of a high growth rate instability. This instability operates in a region 
where magnetic energy density exceeds the rotational energy density. 
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I. INTRODUCTION 



Disk accretion is one of the most fundamental pro- 
cesses occurring in a variety of astrophysical objects like 
accretion powered x-ray pulsars, protostars, cataclysmic 
variables etc. It is observed that the molecular viscosity 
is too small to account for the observed time scales and 
luminosities 1]. To circumvent this problem, Shakura 
and Sunyaev [2[ made a normalization of the coefficient 
of viscosity in the theory of turbulent viscosity proposed 
by Heisenberg Q and introduced a dimensionless param- 
eter ac s H, where c s and H were the sound speed and the 
disk scale height respectively. Though, the a model is 
quite successful in explaining many astrophysical obser- 
vations, a degree of ad-hocness enters the model through 
the parameter a. The presence of Rayleigh criterion of 
hydrodynamic stability, i.e. a flow with specific angular 
momentum increasing monotonically, being satisfied in 
accretion disk, rules out the hydrodynamic origin of tur- 
bulence, though there is literature available [J, l5| in its 
favor. The origin of turbulence is still a question. Dif- 
ferential rotation can lead to an instability provided the 
specific angular momentum decreases outward, a condi- 
tion generally not satisfied in the accretion disk, though 
there are reports suggestingsuch regions being realized 
locally . A recent work Q numerically demonstrated 
the presence of such regions. Though differential rotation 
alone can not lead to instability, inclusion of magnetic 
field alters the scenario altogether. It was well known 
from the classical works of Velikhov 10 and Chandra 
[TT| that a differentially rotating flow with a negative an- 
gular velocity gradient in a weak magnetic field is unsta- 
ble. The presence of this instability in accretion disk was 
established by Balbus & Hawley [12j after three decades 
that the subthermal fields in the presence of differential 
rotation can cause magnetorotational Instability (MRI) , 
which can provide the physical basis for the transport of 
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angular momentum in accretion disk. It has been com- 
putationally shown [l3|, [l4j that MRI is capable of intro- 
ducing turbulence in the accretion disk. It is to be noted 
that MRI is essentially an instability of incompressible 
flow and its growth starts suppressing as the compress- 
ibility seeps in. The earliest attempt to understand the 
effect of compressibility on the growth of MRI was made 
by Blaes & Balbus [15(. They concluded that the 
does not affect the instability. Later on, the behavior of 
compressible MRI in the MHD flows was addressed by 
Kim & Ostriker [l6j]. It was demonstrated that when 
the magnetic field strength is superthermal, the inclu- 
sion of toroidal fields tends to suppress the growth of 
the MRI, and that for quasi-toroidal field configurations 
no axisymmetric MRI takes place in the limit c s — > 0. 
In an another work, Pessah & Psaltis [l?} studied the 
role of toroidal fields in compressible flows. It was again 
demonstrated that the growth rate of MRI is affected in 
the compressible flow and MRI stabilizes in superther- 
mal fields. Clearly, the compressibility plays an impor- 
tant role in the dynamics of instabilities operating in the 
accretion flow. 

Recently, Bonanno and Urpin [l8[ (henceforth paper I) 
studied the effect of compressibility on the instabilities in 
the presence of magnetic field. They considered the mag- 
netic field with non vanishing radial and azimuthal com- 
ponents. A new instability was observed which survived 
for all values of magnetic field, unlike MRI which survives 
only in the weak field limit. The maximum growth rate 
of the reported instability was ~ fi where fi is the rota- 
tion frequency. In a recent work jlfjj , we investigated the 
behavior of above instability for a special case fig < (f^ 
being epicyclic frequency, defined as Vt 2 e = 4fi 2 (l + ^^), 
fi is the rotational frequency) . An instability with a high 
growth rate was observed. In an offshoot of this work, 
we investigate the behavior of instability reported in [lfi | 
in the Keplerian flow. It is to be noted that the only dif- 
ference between the work carried by Bonanno and Urpin 
[lH and us is the inclusion of parameter space which was 
not considered in their study. 

The paper is organized as follows: The section II in- 
vestigates the growth rate of instability. The section III 
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deals with the Results and Discussions. Finally, we con- 
clude by summarizing the new findings. 



II. THE INSTABILITY CRITERIA 



([20j) (see appendix) to find the regions describing in- 
stabilities. This method has been a ppl ied in the field of 
astrophysics very frequently ( flU l21H25j | ) Instability exist 
if any of the conditions written below is satisfied. 



Paper I considered an axisymmetric differentially ro- 
tating system in the presence of a magnetic field. A cylin- 
drical coordinate system (s, tp, z) with s being the radial 
distance from the rotation axis was constructed. Unper- 
turbed system was described by (v r ,v,p,v z ) = (0, sCl, 0). 
Furthermore, il where fl being the angular velocity of the 
astrophysical flow, was taken to be approximately a func- 
tion of s alone; i.e. = fl(s). The isothermal flow for a 
compressible fluid was described by the MHD equations 



i> + (v ■ V)v = -— + g+ t^-(V x B) x B, (1) 
p Anp 



p + V • (pfl) = 0, (2) 



p + v ■ Vp + 7pV • v = 0, (3) 



B - V x (v x B) + 77V x (V x B) = 0, (4) 

V • B = 0. (5) 

where p and v are the density and fluid velocity, respec- 
tively; p is the gas pressure; g is gravity; B is the mag- 
netic field, 77 is the magnetic diffusivity, and 7 is the 
adiabatic index. Magnetic field has non vanishing radial 
and azimuthal components. 

Quasistationary state for a differentially rotating flow 
with magnetic field was considered. Axisymmetric Eu- 
lerian perturbations with space time dependence oc 
exp(crr; — ik ■ r) where k = (fc s ,0, k z ) were introduced 
in the unperturbed accretion disk and the dispersion re- 
lation was obtained for a special case of k.B = after 
neglecting the Ohmic dissipation in the induction equa- 
tion. 

We take the dispersion relation of paper I for determin- 
ing the stability of axisymmetric short wavelength per- 
turbations . The dispersion relation (Eq. (16) of paper 
I) is given by the 



(win) 2 > 0. 



(7) 



er 5 + a 3 {wl + ill) + <r 2 uj 3 Bn + apVL 2 e ujl + pVt 2 e oj 3 



an 
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Equation ([5]) is a polynomial of degree five, so five 
non trivial roots exist. We apply Routh-Hurwitz method 



From the above inequalities, it is clear that the insta- 
bility will exist only when uj%q ^= 

It is to be noted that only those perturbations are con- 
sidered here in which the wavevector is perpendicular to 
the magnetic field, i.e. k ■ B = 0. 

Let us investigate the growth rate of the hydromag- 
netic instability. To calculate the growth rate of this 
instability, it is convenient to introduce the dimensionless 
quantities 



n P 



x 2 Vll ' * 



1 w 



BQ 



c 2 Q3 



x = ks (8) 



The polynomial given by Eq. ((8J) becomes 

T 5 + r 3 (^a; 2 + 1) + T 2 Cx 2 + Tp^x 2 + pQx 2 = (9) 

This equation is solved numerically (see [26| for details) 
by computing the eigen values of the matrix whose char- 
acteristic polynomial is given by equation Q for different 
values of p, £ and £ parameters. 

Since we are interested in the instabilities, only real 
roots (perturbations oc exp(at-ik-r)) are considered. We 
consider the parameter space constrained by £ > £ and 
investigate the growth of instability for different values 
of £ and C, keeping in mind £ > £. It is to be noted that 
this parameter space was not considered by Bonanno and 
Urpin The Fig.(l), (2),& (3) shows the dependence 
of Re (r) for p = 0.3. The value of C is fixed to 1.0, 
2.0 and 3.0 and the values of the parameter £ is varied 
for £ = 0.05, 0.1, 0.2. Similar trend is observed in all the 
cases, i.e. the growth rate of instability increases with 
the decrease in £. For a value of £ = 0.05, growth rate 
of the instability at x 2 — 500 increases from 3.512 e to 
4.6Sl e and 5.50 e for £ = 1.0,2.0,3.0 respectively. As 
x 2 increases further, at x 2 = 700, the growth rate of 
instability increases from 3.8f2 e for ( = 1.0 to 5. 10 e for 
C = 2 to 6.0Sl e for C = 3. It is evident that a vigorous 
instability is operating in the above described parameter 
space. 



III. RESULTS AND DISCUSSION 

We revisit the problem of hydromagnetic instability 
in differentially rotating compressible flows and analyze 
the instability. The parameter space described by ( > £ 
was investigated and diskovered a vigorous instability. 
The behavior of the instability for a typical value of ( 
and £ ( p = 0.3, C = 1.0, 2.0, 3.0, £ = 0.05,0.1,0.2) was 
studied. It is observed that the instability growth rate 
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FIG. 1: The dependence of real parts of the root T on x 2 for 
fi = 0.3, C = 1.0 and £ = 0.2, 0.5, 0.05 



increases with the increase of £ and it also depends on the 
wavenumber. The growth rate of the instability becomes 
very high, as high as (~ 6.012) for x 2 — 700 at £ = 3, 
£ = 0.05). 

In a recent work, Bonanno and Urpin [l8| carried out 
the studies for a compressible flow and showed the pres- 
ence of a new instability. The notable point of their in- 
stability was that it existed for any value of magnetic 
field, the only condition being oj%q ^ 0. The major 
difference between the studies carried out by Kim & Os- 
triker [16j . Pessah & Psaltis [Uj and Bonanno & Urpin 
(l8| is that, in the later case, they considered a field with 
non vanishing radial and azimuthal component. We in- 
vestigated a parameter space and noticed a high growth 
rate instability in that region. Recall that the MRI does 
not exist for k ■ B = 0. Since present study was carried 
out for k ■ B — 0, MRI does not operate in our case. It 
should also be noted that for the incompressible case, i.e. 
c s — > oo, polynomial ((HJ becomes 
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FIG. 2: The dependence of real parts of the root T on x 2 for 
fi = 0.3, C = 2.0 and £ = 0.2, 0.5, 0.05 
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FIG. 3: The dependence of real parts of the root T on x 2 for 
fj, = 0.3, C = 3.0 and £ = 0.2, 0.5, 0.05 



[a 2 + f*%) = 



(10) 



It is clear that the present instability appears only for 
compressible fluid. It is to be noted that the radial com- 
ponent of magnetic field is non zero in the present case. 
Radial field in the presence of differential rotation shears 
into an azimuthal field resulting in an increase in the 
magnitude of Bq. As an example, Desch [27j describes 
that if B^ = B r in a Keplerian disk at the beginning 
of one orbit, by the end of the orbit B^ increases to 
~ 10-B r . After two orbits, B^ increases even more, to 
~ 20B r . This problem, i.e. the increase in the azimuthal 
field in the present case was avoided by treating the state 
as a quasistationary state. When 7/ is small, one can ob- 
tain from Eq. (4) that the azimuthal field grows approx- 
imately linearly with time, 



(11) 



where fl' = dfl/ds, and B v (0) is the azimuthal field at 
t = 0. As long as the second term on the r.h.s. is small 
compared to the first one, and 



t<<T - = s n 



1 B v (0) 



(12) 



stretching of the azimuthal field does not affect signifi- 
cantly the basic state; t v is the characteristic timescale 
of generation of B v . As a result, the basic state can be 
treated as quasi-stationary during the time t <C t v . In 
recent times, Bonanno and Urpin has shown that this 
instability exists in many cases of interest [2l|, 28]. The 
presence of this instability has also been demonstrated 
in the case of protostellar disks [29| also. It is clear that 
the instability discussed in the present report appears in 
many astrophysical situations including protoplanatery 
disk (albeit with small growth rate). 
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A. Difference between MRI and the present 
instability 

This instability is different and independent of the 
magnetorotational instability. The necessary and suffi- 
cient condition for this instability to occur is w| fi ^ 0, 
i.e. B v B s sVL' ^ 0, ft' = therefore, 

1) The instability exists for non vanishing radial and 
azimuthal component of magnetic field whereas the 
poloidal component of magnetic field is important for 
MRI to occur. MRI can exits for vanishing radial and 
azimuthal fields. 

2) The instability exists for any differential rotation, 

i.e. for any sign of 4p whereas MRI exists only when 

the angular velocity in the disk decreases outwards, i.e. 
dn < o 

as 

3) The instability exists even in sufficiently strong mag- 
netic field which suppresses the magnetorotational insta- 
bility. MRI exists only in the weak field limit. 

A special case was considered in the present report 
where wavevector is perpendicular to the magnetic field, 
i.e. k.B = 0. Since the growth rate of MRI is directly 
proportional to the k.B, MRI does not exists in this case. 
This is one of the most important differences between the 
present instability and the MRI. 

B. Where the instability operates 

It is observed that a high growth rate is reached for 
large values of ( > £ (( > 1). If — ft (typical in kep- 
lerian disk), the magnetic energy density is greater than 
the rotational energy density. Therefore, the present in- 
stability will be operating in a regime where "magnetic 
energy density more that the rotational energy density 
" can be realized. For an accretion disk, the rotational 
energy (essentially by definition) is larger than all other 
energies thermal, magnetic, radiation. So, the rotation- 
ally supported accretion disk are out of question to en- 
visage the regions where magnetic energy density is more 
than rotational energy density. One possible realization 
of such regions might appear when one inches towards 
the inner most region of accretion disk, where the disk 
starts truncating/evaporating and the rotational support 
starts getting diminishing and magnetic field takes over. 
A standard pulsar (or even more a magnetar) meets this 
condition. Consider a 10 12 G surface field on a pulsar. 
Several scale heights away from the star the density is 
very small, so the rotational energy density is much less 
than the magnetic energy density. If we consider the low- 
density regions far above the disk then any magnetic bub- 
ble that floats away after reconnection could temporarily 
have a larger magnetic energy density than rotational 
energy density. 



The scenario closest to our regime was considered by 
Bcgclman and Pringlc 30] where magnetic pressure ex- 
ceeds the combined gas + radiation pressure in the disk 
but not the rotational energy density. For magnetohy- 
drodynamic instabilities in accretion disk, the balance is 
achieved via the joint action of rotation, pressure gradi- 
ents and magnetic tension. For the regions mentioned 
above, one can not count on rotation. It is important 
to identify a robust equilibrium about which to perturb, 
perhaps gradients in pressure or radiation pressure could 
help in the inner disk regions. One needs to solve such 
scenario theoretically. We will consider this scenario in a 
future work. 



IV. CONCLUSION 

In the present report, we have demonstrated the pres- 
ence of very high growth rate instability in a differentially 
rotating compressible flows. It has been shown that the 
growth rate of this instability is very high and for a spe- 
cial case considered above, the growth rate increases to 
as high as ~ 6ft for x 2 — 700. It is observed that the 
present instability might be operating in regions where 
magnetic energy density is more than the rotational en- 
ergy density. 
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Appendix A: Hurwitz method 

Let us consider an 5 th order polynomial P(x) 

P(x) = a$x 5 + a^x A + a^x 3 + a 2 x 2 + + oq (Al) 

The Hurwitz theorem states that the above polynomial 
will be unstable if any of the following inequalities is sat- 
ified [13. 

a < , 

Ai = 0,40,3 — 02 < , 
A 2 = 02(0403 - a 2 ) - 04(0404 - a-o) < , 
A 3 = (a 4 ai - a )[a 2 (a4a 3 - a 2 ) - 04(0404 - a )]- 
-a (a 4 a 3 - a 2 ) 2 < , (A2) 
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